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In Ref. Q, we have studied the quantum phase tran- 
sition of the two-dimensional diluted Heisenberg antifer- 
romagnet by using the quantum Monte Carlo method. 
The two main conclusions in the article are i) the critical 
concentration of magnetic sites is equal to the classical 
percolation threshold (p c \ — 0.592746), and ii) the value 
of the critical exponents, such as (3 and iff, significantly 
depends on the spin size S. In gj, Sandvik has presented 
the result of his quantum Monte Carlo simulation for the 
same system only in the case of S = \, in which the 
temperature dependence of the static structure factor is 
studied quite precisely. Based on his data, he claimed 
that the data in Ref. jjj should be seriously affected by 
the finiteness of temperature in our simulation. In addi- 
tion, he calculated the cluster magnetization m c \(L) just 
at p = pd, and claimed that it could remain finite even in 
the thermodynamic limit in the same way as the classical 
(S = oo) case, although quite large finite-size corrections 
to its asymptotic behavior should be considered in the 
quantum case. As we will discuss below, however, his 
data do not conflict with our data nor with our conclu- 
sions. In particular, the data presented in Fig. 1(b) in 
H can be explained from our view point assuming non- 
magnetic ground states. 

As for his first claim on the effect of the temperature, 
we agree that there may be a small underestimation for 
the zero-temperature values of the static structure factor. 
However, even if we replace our data point in Fig. 3 in 
for L = 48 by what we can expect from Fig. 1(a) in 
it would result in only a minor change in the estimate of 
iff. Therefore, the error caused by this underestimation, 
if any, does not seem to be large enough to invalidate our 
conclusions. 

In addition, in Ref. |5j], we have presented the result of 
another finite-size scaling analysis (Fig. 4), in which the 
data at higher temperatures are also used together with 
the one at the lowest temperature in the simulation. In 
this analysis, we have explicitly taken into account the 
finite-temperature effect on the static structure factor by 
introducing the dynamical exponent z. For S = 1, the 
value of the exponent iff coincides with that obtained by 
the former analysis. On the other hand, for S = h, it be- 
comes slightly larger (about 9%) than the other estimate 
based on Fig. 3 in This small difference may be due 
to the finite- temperature effect pointed out in Q , though 
we have not commented on it in [jjj. However, even the 
estimate based on the finite-temperature analysis signif- 


slope = 0.53 


FIG. 1. Log-log plot of cluster magnetization (Ref. M). 
The solid line is obtained by least-squares fitting for L > 28. 

icantly differs from those for S > 1. Therefore, we think 
our second conclusion on the non-universal behavior can 
not be excluded. 

It was also suggested (|], based on the calculation of 
the magnetization of the largest connected cluster for 
each sample, that all the critical exponents, including 
iff and j3, are equal to the classical values. The data were 
plotted against x = 1/L D / 2 , and appeared to suggest a 
finite value in the thermodynamic limit. From our point 
of view, however, his data can be interpreted as follows. 
The quantity m 2 j(L) defined in Ref. ||] should be asymp- 
totically scaled as 


S s (L,0,Pd)/L 


2D-d 


(1) 


where S s (L,T,p) is the static structure factor (Eq. (4) 
in Ref. jl]]) and L D is the average number of spins in the 
largest connected clusters in a system of linear size L. Ac- 
cording to our results pj], on the other hand, S S (L, 0,p c \) 
diverges as L 2D ~ d ~ a with a — 0.52 for S — |, and con- 
sequently m^L) should be scaled as L~ 52 . In Fig. [I], 
we show a log-log plot of the cluster magnetization data 
by Sandvik B . One sees that the data for L > 28 clearly 
scales as L~^° 3 , which is fully consistent with our results 
mentioned above. 

The authors are grateful to A. W. Sandvik for allowing 
us to use his numerical data. 
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